Suppresion of Bloch oscillations by weak disorder in semiconductor superlattices 
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We investigate the dephasing dynamics of Bloch oscillations in semiconductor superlattices by 
means of a very simple model including disorder and applied electric fields. A thorough numerical 
study of our model allows us to claim that small, unintentional well width fluctuations can be 
responsible for fast dephasing of Bloch oscillations at low temperatures. We show that the lifetime 
of Bloch oscillations is controlled by a characteristic time which depends on the degree of disorder 
and is independent of the electric field. This result is further supported by the excellent agreement 
between our model calculations and several recent experiments, and leads to specific new predictions. 
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Dynamical effects in quantum- well semiconductors su- 
perlattices (SL's) are the basis for designing ultra-high 
speed electronic devices, as have been recently proposed 
0. This idea of semiconductor SL's operating at tera- 
hertz frequencies was already suggested a long time ago 
by Esaki and Tsu ||], who argued that electrons should 
undergo periodic Bloch oscillations (BO's) Under 
an applied electric field F, provided that the interband 
coupling is negligible, electrons oscillate in real space as 
well as in k space with a characteristic period given by 
TB — 2TTfi/eFd, d being the spatial period of the SL 
The amplitude of BO's in real space is A ~ A/2eF, where 
A is the minibandwidth. The coherent carrier motion 
is thus restricted to a region of length 2 A. This peri- 
odic motion persists until the Bloch electron loses energy 
gained from the field through scattering processes. Re- 
ports of unambiguous experimental evidences for BO's in 
GaAs-Gai_2:Al2;As are presently available [^-10|. 

Inelastic scattering by phonons, deviations from SL's 
perfect periodicity due to unintentional imperfections, in- 
traband scattering, interminiband transitions, and scat- 
tering by impurities severely reduce the quantum coher- 
ence required for the observation of BO's. Indeed, the 
scattering time r must be larger than the Bloch period tb 
and therefore the electric field must be larger than certain 
critical electric field Fc Q. However, even in the most 
favorable experimental conditions t is not much larger 
than Tb and thus only a few BO's are usually observed. 
The origin of such loss of quantum coherence in actual 
devices is far from understood and, at present, there is 
much debate about the role played by different scattering 
mechanisms in those processes. In this regard, Plessen et 
al. 1^ found that quantum coherence is lost after few 
BO's in 30 A GaAs/30A Gao.7Alo.3As SL's, which was 
attributed to scattering by LO phonons. On the other 



hand, theoretical studies point out that under most ex- 
perimental conditions interminiband transitions are neg- 
ligible and, consequently, cannot be responsible for the 
signal decay |^]. Furthermore, Plessen et al. Q conclude 
from their experimental results that the critical electric 
field Fc is higher for SL's with A larger than the energy 
of LO phonons, E'lo = 36meV. They explain this depen- 
dence by assuming that LO phonon emission is excluded 
when A < Slo- On the contrary, Leisching et al. jl^] 
detected up to six BO's but they did not observe any 
sign of a phonon threshold in SL's with A ranging from 
13 up to 46 meV. These authors argued that the reduced 
sample quality of Ref. could be the responsible for the 
threshold. 

From the above discussions, it becomes clear that un- 
derstanding the interplay between the electric field and 
the imperfections of the SL's is crucial to elucidate the 
discrepancies among different groups, either to pinpoint 
its relevance or to exclude it. As far as we know, however, 
a complete study of the effects of interface roughness on 
the Bloch oscillations dynamics is currently lacking. In 
this letter, we introduce a theoretical model for imper- 
fect GaAs-Gai_2;Alj;As SL's that successfully accounts 
for the experimental results. We study the dynamical 
behavior of these disordered SL's subject to a dc elec- 
tric field by measuring the position of the centroid of the 
wavepacket and by means of the time dependent inverse 
participation ratio (IPR), to be defined below. These 
quantities will allow us to conclude that the assumption 
of weak disorder is enough to explain all the available 
experimental data, thus firmly connecting the dephasing 
of BO's to the quality of the sample. 

Interface roughness appearing during growth in actual 
SL's depend critically on the growth conditions For 
instance, protrusions of one semiconductor into the other 
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cause in-plane disorder and break translational invari- 
ance parallel to the layers. To describe local excess or 
defect of monolayers, we allow the quantum well widths 
to fluctuate uniformly around the nominal values; this 
can be seen as substituting the nominal sharp width by 
an average along the parallel plane of the interface im- 
perfections. Our approximation is valid whenever the 
mean-free-path of electrons is much smaller than the in- 
plane average size of protrusions as electrons only see 
micro-quantum- wells with small area and uniform thick- 
ness . In each micro-quantum- well presents a slightly 
different value of its thickness and, as a consequence, 
resonant coupling between electronic states of neighbor- 
ing GaAs layers is decreased. Therefore, in the follow- 
ing we will take the width of the nth quantum well as 
a(l -I- Wen), where 14^ is a positive parameter measuring 
the maximum fluctuation, e„'s are distributed accord- 
ing to a uniform probability distribution, P(e„) = 1 if 
|e„| < 1/2 and zero otherwise, a is the nominal quantum 
well width. For clarity we assume that the barrier width 
b is the same in the whole SL, although we have checked 
that this assumption can be dropped without changing 
our conclusions. 

For our present purposes, it is enough to focus on elec- 
tron states close to the conduction-band edge and use the 
effective-mass approximation. The envelope-functions 
for the electron wavepacket satisfies the following quan- 
tum evolution equation 
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where x is the coordinate in the growth direction and V 
is the SL potential at flatband. We have considered a 
constant effective-mass m* at the F valley for simplicity, 
but our numerical results should qualitatively describe 
actual SL's with position-dependent effective-mass. 

We study the quantum dynamics of an initial Gaussian 
wavepacket 
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where the mean kinetic energy is {E) — h k'^/2m* and a 
measures the width of the electron wavepacket. We stress 
that, according to Bouchard and Luban the dynami- 
cal behavior of this initial state is similar to that of more 
realistic functions. The solution of Eq. (Q) is accurately 
obtained using the Cayley's form for the finite difference 
representation of its formal solution Once the 

solution is obtained, we evaluate the position of the cen- 
troid of the wavepacket as 



/oo 
dx {x- XQ)\^{x,t)[-' 
-OC 



(3) 



which should display BO's. Moreover, we use the time- 
dependent inverse participation ratio (IPR), defined as 
the second moment of the probability density 
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to evaluate the spatial extent and the degree of localiza- 
tion of electronic wavepackets. We note that delocalized 
states present small IPR (in the ballistic limit, without 
applied field, it vanishes with time ast~^), while localized 
states have larger IPR. 

We have considered the same parameters as those of 
the SL's used in previous experiments |^,|o|. In particu- 
lar, we present here results for the first one of these SL's, 
i.e., 100 periods of 30 A GaAs and 30 A Gao.7Alo.3As 
Q . Samples are labeled according to their period length 
d — a + b, namely 6OA SL. Similar results are obtained 
with the other SL's like the 84, 97, or 128 A (6 = 17 A, 
a = 67, 80 and 111 A, respectively), i.e., the ones re- 
ported by Leisching et al. |l^] , although we do not present 
here these results for brevity. We have straightforwardly 
calculated the miniband-width for the 6OA SL obtaining 
A = 90meV, being larger than Ej^q. We study apphed 
electric field in the ranges from 5 up to 20kV/cm. The 
fluctuation parameter runs from W = (perfect SL's) 
up to VF = 0.20 (strongly disordered SL's). This max- 
imum value considered here represents excess or defect 
of a few monolayers. This value is above the degree of 
perfection now achievable with MBE, so that realistic re- 
sults are comprised within this range and we do not need 
to analyze stronger disorder values. 

Figure |l| displays the centroid position of the 
wavepacket in the 60 A SL for F = lOkV/cm and differ- 
ent values of the unintentional disorder. The initial Gaus- 
sian wavepacket is located in the centermost quantum 
well with a = 300 A and fco = . In Fig. |l|(a) , for an or- 
dered SL, we observe the occurrence of very well defined 
BO's with amplitude 2A ~ 900 A and period tb = 0.7 ps, 
in excellent agreement with the theoretical predictions 
2A — 900 A and tb — 0.69 ps. Notice, however, that the 
perfect oscillatory pattern detected in periodic SL's (up- 
per panel) is progressively destroyed upon increasing the 
degree of disorder as seen in the rest of panels in Fig. [|, 
for W = 0.01, 0.03, 0.05, 0.10, and 0.20. It is most impor- 
tant to mention here that the results do not depend on 
the particular realization of disorder. We note that those 
values correspond, if we assume that a monolayer width 
of this type of SL's is about 3A, to a maximum excess of 
defect of less than one monolayer {W = 0.01 and 0.03), 
one monolayer {W — 0.05), two monolayers {W — 0.10) 
and four monolayers {W — 0.20). The disorder induces a 
decrease of the amplitude of the oscillations and, besides, 
it produces a progressive dephasing comparing with the 
ideal perfect case [Fig 0(a)] . In the strong disorder case 
no signs of BO's are found. This fact can be explained by 
the absence of translational invariance at flatband and, 
consequently, by scattering of electrons with the random 
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potential. Similar results are obtained with the SL's re- 
ported by Leisching et al. [0. We have to note in this 
respect that for those samples we have observed a slower 
decay of BO's, the slowest evolution corresponding to the 
SL with the narrowest miniband. This can be understood 
by noting that the amplitude of BO's is proportional to 
the miniband width: Thus, the coherent motion of elec- 
trons takes place in a smaller region when decreasing the 
mini-band-width and, therefore, carriers are less influ- 
enced by disorder. This result is in very good agreement 
with the experimental observations pO| . 
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FIG. 1. Centroid of an initial Gaussian wavepacket with 
fco = and o- = 300 A as a function of time in 60 A SL's. 
The applied electric field is F = lOkV/cm. From top to bot- 
tom (a) W = 0, (b) 0.01, (c) 0.03, (d) 0.05, (e) 0.10, and 
(f) 0.20. The values of the amplitude 2A = 900 A and Bloch 
period r = 0.7 ps., for the perfect SL's (a), are in excellent 
agreement with the theoretical predictions. 

We can achieve a better resolution of the BO's period 
and the influence of the disorder by means of the IPR. 
The upper panel of Fig. || presents the results for the IPR 
of the 60 A SL when the initial Gaussian wavepacket is 
located in the centermost quantum- well with ct = 20 A 
and fco = 0. The electric field is F = lOkV/cm. In 
the absence of imperfections, the IPR displays a periodic 
pattern with marked peaks at times i„ = nr^ , where n is 
any arbitrary, nonnegative integer and tb = 0.7 ps. This 
means that the initial localized state is recovered after 



this time. It is most important, to assess the accuracy 
of our calculation, to mention that the numerical value 
of the IPR at maxima is exactly the same than that ob- 
tained from (|) and (|), that is, IPR(O) = l/(20Fcr) = 
0.01 A~^. Results corresponding to disordered SL's with 
the same initial conditions as before are shown in the 
remaining panels of Fig. ^ confirming that BO's pro- 
gressively disappear on increasing the degree of disorder. 
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FIG. 2. IPR vs time for an initial Gaussian wavepacket 
with fco = and a — 20 A, subject to an electric field 
F = lOkV/cm in 60 A SL's. From top to bottom W = 0, 
0.01, 0.03, 0.05, 0.10, and 0.20. 



From the above results we are led to the conclusion 
that there exists a characteristic scattering time r^is af- 
ter which BO's are destroyed by disorder. Moreover, it 
is readily observed in Fig. || that Tdis decreases upon in- 
creasing the degree of disorder. However, the above re- 
sults have been obtained for a fixed value of the electric 
field, but clearly a meaningful definition of the scattering 
time should be independent of the value of the electric 
field. To check the validity of the introduced Tjis we have 
studied the IPR for different values of the applied electric 
field at a given degree of disorder. Representative results 
are presented in Fig. ||for W — 0.03 (on average less than 
one monolayer) and F = 5, 10, 15 and 20kV/cm. From 
this plot we can estimate that rjis — 2.5 ps for all val- 
ues of the electric field. Thus, this scattering time plays 
the same role as the scattering time arising from inelastic 
interactions, in the sense that tb must be kept smaller 
than Tdis to observe BO's. Interestingly, this value is the 
same as that obtained in the experiments of Plessen et 
al. [|| The scattering time increases when the miniband- 
width decreases, for the same amount of disorder, and 
values obtained with our model turn out to be perfectly 
consistent with all the experimental values ||7|, p^ , p^ . 
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FIG. 3. IPR vs time for an initial Gaussian wavepacket 
with the same parameters as in Fig. 2, placed in a 60 A SL's 
with W = 0.03. From top to bottom F = 5, 10, 15, and 
20kV/cm. 



erating devices demand higher electric fields. Therefore, 
for sufficiently high fields, the region where coherent car- 
rier motion takes place, namely 2 A — A/eF, is compa- 
rable to the SL's period d. In such a situation, the in- 
plane disorder is no longer well described by an ensemble 
of different quantum wells as we have proposed because 
the wavepacket only would see one quantum well in our 
model. Therefore more theoretical work is needed to in- 
vestigate the role of imperfections and other dephasing 
mechanisms like excitonic effects in the design of fu- 
ture, shorter period devices. 

The authors thank with great pleasure collaborative 
and illuminating conversations with Fernando AguUo- 
Rueda, Karl Leo and Gintaras Valusis. This work 
has been supported by CICYT (Spain) under project 
MAT95-0325. 



To conclude, we have been able to firmly connect BO's 
suppression and dephasing in actual SL's to small de- 
viations from exact flatness at well-barrier interfaces. 
Specifically we have shown that an average degree of 
imperfection of less than a monolayer suffices to ex- 
plain quantitatively the experimental results in |p|,p!o[. 
Whereas the initially localized state is recovered after 
time tb in the case of perfect {W = 0) SL's (regular 
behavior), any degree of disorder due to imperfections 
during growth leads to the disappearance of BO's after 
a few oscillations: The higher the degree of disorder the 
faster the vanishing of BO's. The very good agreement 
with previous experiments points out the crucial role of 
imperfections in the dynamics of actual SLs driven by 
electric fields. Most importantly, we have been able to 
define a characteristic scattering time Tdis, independent 
of the electric field, after which BO's cannot be detected, 
this being a specific prediction of our model that can be 
checked in experiments. In other words, for the BO's 
to be observed in actual SLs, the apphed electric field 
must be larger than some critical electric field given by 
eFdisd = 27r7i/Tdis. The existence of such a critical field is 
evidently very important from the viewpoint of practical 
applications of our results, i^dis is directly related to the 
degree of disorder present in the sample and decreases 
upon increasing the quality of the sample, i.e., it is an 
excellent parameter to asses the performance of epitaxial 
growth techniques. 

As we have seen, the main conclusion of this work, 
that the importance of disorder in the transport prop- 
erties of SL's has to be underlined, in contrast with the 
general belief than the high quality of actual SL's al- 
lows to neglect disorder as a second order effect. It has 
to be kept in mind that such an average of disorder of 
less than a monolayer is currently unavoidable, more so 
when preparing such long SL's (100 periods) as we have 
considered. We note, however, that high-frequency op- 
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